Off Resonant Pumping for Transition from Continuous to Discrete Spectrum and 
Quantum Revivals in Systems in Coherent States 
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We show that in parametrically driven systems and, more generally, in systems in coherent states, 
off-resonant pumping can cause a transition from a continuum energy spectrum of the system to 
a discrete one, and result in quantum revivals of the initial state. The mechanism responsible for 
quantum revivals in the present case is different from that in the non-linear wavepacket dynamics of 
systems such as Rydberg atoms. We interpret the reported phenomena as an optical analog of Bloch 
oscillations realized in Fock space and propose a feasible scheme for inducing Bloch oscillations in 
trapped ions. 
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I. INTRODUCTION 

Parametrically driven systems have been studied ex- 
tensively in quantum and nonlinear optics. Such systems 
have resulted in the production of twin beams of pho- 
tons with high degrees of entanglement |^ . These beams 
of light have been used extensively in many applications 
such as in connection with Bell's inequalities and EPR 
paradox , quantum cryptography ^ , teleportation , 
quantum imaging |^ and in general, in connection with 
the subject of quantum information processing The 
standard Hamiltonian describing such a parametric inter- 
action is known to have continuous spectrum which then 
results in the exponential growth of (say) the energies in 
the signal and idler modes. The state of the system at 
any given time is, in fact, the so-called Perelomov coher- 
ent state 0. 

In this communication, we report unusual features in 
the dynamics of systems in coherent states. We show how 
changing the detuning of the field pumping the system 
can cause a transition from a continuum energy spec- 
trum of the system to a discrete one. This transition 
from continuum spectrum to discrete spectrum results in 
quantum revivals of the initial state. 

This study has been motivated primarily by recent 
work on the optical realization of Bloch oscillations ||] 
and Wannier-Stark ladders These have been recently 
observed in ultracold atoms trapped in optical potentials 
[0, and in fields propagating in waveguide arrays with 
a linear bias so that the refractive index of a guide is 
proportional to its spatial location Even the effect 
of nonlinearities in the waveguide on Bloch oscillations 
has been studied. The underlying mathematical equa- 
tion which leads to Bloch oscillations has the structure 



iCn = nAC„ + /3(C„+i + Cn-l) 



(1) 



In the absence of the term nAC„, an initial excitation 
in one of the C's gets delocalized , i.e, it spreads over 
all the other sites. However, if A ^ 0, then the sys- 
tem returns to the original site. The linear dependence 
on the site index in the term C„ is very critical for the 



revival of the initial state. In the waveguide example, 
the linear bias produces the linear n dependence of the 
term nAC„. The optical realization of analogs of Bloch 
oscillations is based on situations where the dynamics 
could be described by equations like (|^). We note that 
many coherently driven systems have dynamics charac- 
terized by equations like (|l|). Thus, remarkably enough, 
systems in coherent states can exhibit analog of Bloch 
oscillations. In the case of the radiation field, the in- 
dex n would refer to the occupation number of a given 
Fock state. Thus the standard measurement of the pho- 
ton number distribution can be used to study the analog 
of Bloch oscillations. Clearly this opens up the possibil- 
ity of studying Bloch oscillations in a very wide class of 
driven systems. It is however important to have A 7^ 0, 
that is that the pumping should be off-resonant. Note 
also that the quantum revivals considered in this paper 
should be distinguished from the well studied quantum 
revivals in the nonlinear wavepacket dynamics of many 
systems such as Rydberg atoms [[l2| and coherent states 
propagating in a Kerr medium p3|T 



II. PARAMETRICALLY DRIVEN SYSTEMS 

Consider a very commonly used system in optical 
physics. The Hamiltonian describing the parametric in- 
teraction of two modes a and b with frequencies uia and 
LOi, respectively, is given by 



H 



(2) 

Here G is the coupling constant which depends on the 
nonlinearity of the medium and the amplitude of the 
pump. It is assumed that the pump frequency ujp — 
uja + LOb. The operators a and b satisfy Bosonic commu- 
tation relations. If the modes a and b are in vacuum state 
at time t = 0, then the state at time t is 



\i^{t)) = exp{-it(Gat6t + G*a6)}|0,0) 



(3) 



which is the so-called Perelomov coherent state or the 
squeezed vacuum state . As is well known, the Heisen- 
berg operator a, for example, evolves as 



1 



i{t) ^ a(0) cosh(|G|i) - ib^O)-^ sinh(|G|t) (4) 



and thus the energy in the mode a continues to grow: 
{a^t)a{t)) = sinh^(|G|t). The probability po(t) of find- 
ing the system in the initial state is 

= |(0,0|exp{-ii(Ga^6t + G*a6)}|0,0)p 

= cosh~2(|G|i)-4e-2|G^I*, as i ^ oo. (5) 

Thus the probability of returning to the initial state de- 
creases monotonically and becomes exponentially small 
in the long-time limit (Fig. 1). Clearly, there is no quan- 
tum revival of the initial state. 

The corresponding probability p„(i) {n ^ 0) of finding 
the system in a state with n photons in each mode, is 
given by 



p„(t) = cosh"2(|G|i) tanh2"(|G|i); n ^ 0. 



(6) 



Note that for n ^ Q, Pn{t) increases with t initially and 
peaks at \G\t = tanh^^ yjn/{n+ 1). The peak value of 
Pn{i) has the expression 



Vn{t)\ 



peak 



(7) 



which depends only on n and decreases monotonically 
with n. 




|G|t 

FIG. 1. Evolution of a parametrically driven system when 
the pumping is resonant (A = 0). The system dynamics 
is governed by the Hamiltonian given by Eq. (2). Here, 
Pnit) is the probability that the system, starting from vac- 
uum, reaches a state at time t with n photons in each mode. 
Plots of Pn{t) as a function of \G\t are shown for n = 0, 1, 
2 and 3. The plots are identified by their peak values which 
decrease as n increases. Note that the probability of return- 
ing to the initial state (n — 0) decreases monotonically with 
time. 



The exponential growth of energies in the signal and 
idler modes and the subsequent absence of quantum re- 
vival for the initial state are connected with the fact that 
the Hamiltonian admits a continuum of eigenvalues. 
The corresponding eigenstates are delta function normal- 
izable and signify scattering states fl^ : 



H\i,ix)) = mx)) 

(VXA')IV^(A)) = Six - A'); A, A' real. 
IV'(A)) = ^d„|n,n), 



(8) 



where the coefficients dn are expressed as hypergeometric 
functions of the form 2^1 (a, b; c, z). 

Writing ]?/'(<)) = X^n ^^'^ substituting in 

the Schrodinger equation ih\^p) = H\^) with H as given 
in (|), we get 



iCn = GnCn-i +G*{n + l)G„+i. 



(9) 



Thus in terms of Fock states we have nearest-neighbor 
coupling with the coupling itself depending on the index 
n. The solutions (^)and imply that with increasing 
time, an initial excitation in the n = state spreads over 
all the states, or in the language of sohd state physics, 
the excitation delocalizes. 

We now consider an important physical situation 
where the Hamiltonian (0) is modified to 



H 



— {a^a- 



b^) + h{Ga''b'' + G*ab). 



(10) 



This will be the case if the pump field is detuned, i.e, 
ujp LUa + ^b- The amplitudes will satisfy Eq. (|^) with 
an extra term nAG„ on the right. The state at time t can 
be obtained by using the general disentangling theorem 
||15| for the SU(1,1) operators K+ ^ a'^b^ K_ = ab and 
K3 ^ {a^a + b^+l)/2: 

exp{it{GK+ + G*K_ + AK3)} 
= exp{itGT+K+) exp{(lnr3)Js:3} exp{itG*T_K. 

where 

/ tA ' 

Fs = ( cosh/? — z— sinh/3 



(11) 



F+ = 



2sinh/3 



2(3 cosh/3-iiAsinh/3' 

/3o = VAV4-|G|2. 

The probabilities are now given by po(t) 
Pn{t) - |F3||F+|2" where 



(12) 
and 



IF., 



|F+P = 



-i\G\/(3ofsm^Pot 

sin^ (3ot 
sin^/3oi+(/3o/|G|)2 



(13) 
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if the detuning A and the couphng constants are such 
that /3o (defined by Eq. (|l^)) is real. Clearly the sur- 
vival probability of the initial state is now an oscillatory 
function of time t and po{t) is unity whenever Pot is an 
integer multiple of tt, i.e. 

(A^/4 - |Gp)i^ = n^TT^, where n is an integer. (14) 

That is, for these times, the system returns to the origi- 
nal state |0, 0) up to an overall phase factor and there is 
a quantum revival (Fig. 2). 



where q = a^a — b^b is an integer constant of mo- 
tion. Thus the periodic behavior is connected with a 
transition in the nature of energy levels, i.e, the system 
has gone from a continuous spectrum to a discrete spec- 
trum whose energy levels are equidistant. Note that 
the level separation is 2h(3o which further validates our 
earlier observation that the revival period is given by 
7r/^A2/4 - |Gp (see Eq. (|lj)). Comparison with the 
standard Bloch problem suggests that A plays the role of 
the electric field whereas the magnitude of G determines 
the coupling to the nearest neighbors. 




FIG. 2. As in Fig. 1 except that the pumping is 
off-resonant ( A 7^ 0) . The corresponding Hamiltonian is given 
by Eq. (0). Plots of p„(f) are shown for (a) A/2lG| = 1.01 
and (b) A/2|G| = 1.10. 



We thus find the important result that our parametric 
system exhibits analog of Bloch oscillation under the con- 
dition of non-zero detuning such that A > 2|G| and that 
the survival probability becomes unity for times defined 
byEq. (0). 

It is interesting to note that the peak values of pn (t) 
(n 7^ 0) are still given by Eq. (Q) although their loca- 
tions, at /3o< = sin^^(^/n/3o/|G|), will now depend on the 
detuning A. 

Under the condition A > 2|G|, the Hamiltonian ( p^ ) 
admits a discrete set of eigenvalues given by 

= h^oiq + 2m + 1) - hA/2; m = 0, 1, 2, ... (15) 



III. OSCILLATOR DRIVEN BY A DETUNED 
COHERENT FIELD 



All the above has analogs in the theory of coherent 
states for a single mode of the radiation field. In fact, 
coherent states will exhibit some new interesting physics 
if the detuning is nonzero. Consider the behavior of an 
oscillator driven by a detuned coherent field: 



H = h{Aa^a 



a). 



(16) 



If the system is in a vacuum state at t — 0, then the state 
at time i is a coherent state: 



\ip{t)) = exp 
'1 - 



-iAt 



sin At) 



iA 



(17) 



Consider first the case of zero detuning, A = 0. In this 
case, the amplitude —iet of the coherent state continues 
to increase as t increases. The peak of the photon num- 
ber distribution continues to move away from zero as t 
increases and the probability of being in the initial state 
decreases continuously and becomes exponentially small: 
Po(i) = exp(-|e|2i2). 

For non-zero detuning, on the other hand, one obtains 



p„it) = \{n\m)\ 



1 \a 



2n 



n > 0; 
\a\ = \e\t 



|sin(At/2)| 



Ai/2 



(18) 



Note that the peak of Pn{t) occurs at n ~ |ap = 
(2e/A)2sin2(At/2) which is an oscillatory function of 
time and vanishes whenever At is an integer multiple of 
27r. At such times, Pn{i) = Sn,o, i-e, the system returns 
to the initial state |0) (Fig. 3) up to a phase factor. Note 
that the period does not depend on the strength of the 
driving field. 
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FIG. 3. Evolution of an oscillator driven by a coherent field 
of amplitude e and detuning A. The corresponding Hamilto- 
nian is given by Eq. (^). Plots of Pn{t) as a function of 
\e\t are shown for ti = 0, 1 and 2 when (a) A/2|e| = 0, (b) 
A/2ie| = 0.5, (c) A/2|e| = 1.0 and (d) A/2|e| = 2.0. The plots 
are identified by their peak values e~"n"/n!, which decrease 
as n increases. All other notations are as in Fig. 1. 



We also note that for A ^ 0, the eigenvalues of the 
Hamiltonian (p^) are equally spaced discrete states with 



The spacing is independent of e which leads to a periodic 
behavior independent of e. On the other hand, if A = 0, 
then H — Ti{ea'^ + e*a) has a continuum of eigenvalues 
with the eigenfunctions satisfying the usual delta func- 
tion condition {'ijj{x)\'ip{x')) = S{x — x'). Thus here again 
we see the transition from continuum eigenvalues to the 
discrete set and this transition is responsible for the 
revival of the initial state. 



IV. QUANTUM REVIVALS IN ION TRAPS 

We now show that this simple behavior can be easily 
probed, as for example, for a trapped ion like the one 
used by Wineland and co-workers [Q. Consider an ion 
in a deep trap cooled to its ground state. Let it interact 
with two Raman fields Ei and E2 with frequencies Wi 
and 072 such that wi — 0^2 — 7^ (Fig. 4). 
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FIG. 4. Raman-like off-resonant transitions among the vi- 
brational states of an ion cooled to its electronic ground state 
in a trap. The ion is acted upon by two fields with ampli- 
tudes El and E2 and frequencies uji and 0^2 respectively. The 
separation between the relevant vibrational levels is denoted 
by v so that the detuning A = uii — UJ2 ~ i^- The ion would 
return to its ground state whenever At = 2nn, where n is an 
integer. 



In such a case the effective Hamiltonian describing the 
ionic motion will be of the form iJii^jO exp{i(a;i — lu2 — 
i')t} + h.c. This Hamiltonian is equivalent to (|l6|). Here 
a is the annihilation operator for the vibrational motion 
of the ion. A measurement of the vibrational distribu- 
tion will exhibit the analog of Bloch oscillation provided 
wi — a;2 — 7^ and the ion would return to its ground 
state whenever (uJi — 102 — i^)t = 2nn, where n is an inte- 
ger. 
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V. CONCLUSION 

In conclusion, we have introduced ofiF-resonant pump- 
ing of coherent systems as a new mechanism for the opti- 
cal analog of Bloch oscillations. In our scheme, such os- 
cillations occur in Fock space. We showed that for para- 
mctrically driven systems and systems in coherent states, 
the detuning of the pump changes the energy spectrum of 
the system from continuum to discrete , and revives the 
initial state. Finally, wc proposed a scheme for inducing 
Bloch oscillations in trapped ions. 
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